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^a_>^ , Abstract. Two cornerstones of the Kasteleyn method are: 1. rewriting the Ising partition function as 

the dimer partition function, that is, the generating function of the perfect matchings, and 2. expressing 
' the dimer partition function of planar graphs as the determinant. This paper initiates the 3-dimensional 

^ ^ Kasteleyn method. 

, We show that the weight enumerator of any binary linear code is polynomial reducible to the permanent 

' of a 3-dimensional matrix (3-matrix). In analogy to the standard (2-dimensional) matrices we say that a 

3-matrix A is Kasteleyn if signs of its entries may be changed so that, denoting by A' the resulting 3-matrix, 
we have per(j4) = det(A'). We show that in contrast with the 2-dimensional case the class of Kasteleyn 
3-matrices is rich; namely, for each matrix M there is a Kasteleyn 3-matrix A so that per(M) = per(j4) = 
det(yl'). 

' In particular, the dimer partition function of a finite 3-dimensional cubic lattice may be written as the 

• determinant of the vertex-adjacency 3-matrix of a 2-dimonsional simplicial complex which preserves the 

, natural embedding of the cubic lattice. 

^: 

1. Introduction 

^ , The Kasteleyn method is a way how to calculate the Ising partition function on a finite graph G. It goes as 

CNJ ' follows. We first realize that the Ising partition function is equivalent to a multivariable weight enumerator 

I of the cut space of G. We modify G to graph G" so that this weight enumerator is equal to the generating 

function of the perfect matchings of G', perhaps better known as the dimer partition function on G". Such 
generating functions are hard to calculate. In particular, if G" is bipartite then the generating function of 
CSj ' the perfect matchings of G' is equal to the permanent of the biadjacency matrix of G' . If however this 

permanent may be turned into the determinant of a modified matrix then the calculation can be successfully 
carried over since the determinants may be calculated efficiently. Already in 1913 Polya [6 asked for which 
non-negative matrix M we can change signs of its entries so that, denoting by M' the resulting matrix, 
we have per(M) = det(M'). We call these matrices Kasteleyn after the physicist Kasteleyn who invented 
I the Kasteleyn method [21 El HJ [10]. Kasteleyn [4] proved in 1960's that all biadjacency matrices of the 

H ' planar bipartite graphs are Kasteleyn. We say that a bipartite graph is Pfaffian if its biadjacency matrix is 

- - -' Kasteleyn. The problem to characterize the Kasteleyn matrices (or equivalently Pfaffian bipartite graphs) 

was open until 1993, when Robertson, Seymour and Thomas |7] found a polynomial recognition method and 
a structural description of the Kasteleyn matrices. They showed that the class of the Kasteleyn matrices 
is rather restricted and extends only moderately beyond the biadjacency matrices of the planar bipartite 
graphs. 

In this paper we carry out the Kasteleyn method for general binary linear codes. We show that the weight 
enumerator of any binary linear code is polynomial reducible to the permanent of the triadjacency 3-matrix 
of a 2-dimensional simplicial complex. In analogy to the standard (2-dimensional) matrices we say that a 
3-dimensional non- negative matrix A is Kasteleyn if signs of its entries may be changed so that, denoting 
by A' the resulting 3-dimensional matrix, we have per(^) = det{A'). We show that in contrast with the 
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2-dimcnsional case the class of Kastcilcyn 3-dimensional matrices is rich; namely, for each 2-dimcnsional 
non- negative matrix M there is a 3-dimensional non- negative Kasteleyn matrix A so that per(M) = per(^). 
Finally we conclude with some remarks directed towards possible application for the 3-dimensional Ising and 
dimer problems. 

1.1. Basic definitions. We start with basic definitions. A linear code C of length n and dimension d over 
a field F is a linear subspace with dimension d of the vector space F". Each vector in C is called a codeword. 
The weight w{c) of a codeword c is the number of non-zero entries of c. The weight enumerator of a finite 
code C is defined according to the formula 

Wc{x) :=^a;'"("^ 

A simplex X is the convex hull of an affine independent set V in R'^. The dimension of X is |T^| — 1, 
denoted by dimX. The convex hull of any non-empty subset of V that defines a simplex is called a face 
of the simplex. A simplicial complex A is a set of simplices fulfilling the following conditions: Every face 
of a simplex from A belongs to A and the intersection of every two simplices of A is a face of both. The 
dimension of A is max {dim X|X e A}. Let A be a d-dimensional simplicial complex. We define the incidence 
matrix A = (Aij) as follows: The rows are indexed by {d — 1) -dimensional simplices and the columns by 
d-dimensional simplices. We set 



Aij :— 



1 if (d — l)-simplex i belongs to d-simplex j, 
otherwise. 



This paper studies 2-dimensional simplicial complexes where each max:imal simplex is a triangle or an edge. 

Wc call them triangular configurations. We denote the set of vertices of A by V{A), the set of edges by 
E{A) and the set of triangles by r(A). The cycle space of A over a field F, denoted ker A, is the kernel of 
the incidence matrix A of A over F, that is {xj^a; = 0}. 

Let A be a triangular configuration. A matching of A is a subconfiguration M of A such that ti D t2 
does not contain an edge for every distinct ti,t2 & T{M). Let A be a triangular configuration. Let M be a 
matching of A. Then the defect of M is the set E{T)\E{M). The perfect matching of A is a matching with 
empty defect. We denote the set of all perfect matchings of A by 'P(A). Let w : T(A) R be weights of the 
triangles of A. The generating function of perfect matchings in A is defined to be Pa{x) = X]p£p(A) x^^^\ 
where w{P) := I]t£pw(i)- 

A triangular configuration A is tripartite if the edges of A can be divided into three disjoint sets Ei, E2, E3 
such that every triangle of A contains edges from all sets £^2, -E3. We call the sets Ei,E2,Ez tripartition 
of A. 

The triadjacency 3-matrix A{x) = {aijk) of a tripartite triangular configuration A with tripartition 
E\,E2,E^ is the \Ei\ x \E2\ x \E^\ three dimensional array of numbers, defined as follows: We set 



^w{t) j£ g. ^ Ei,ej e E2, Cfc e £^3 form a triangle t with weight w{t), 
otherwise. 



The permanent of a n x n x n 3-matrix A is defined to be 



per(A)= \{aia^{i)a2{i)- 
The determinant of a n x n x n 3-matrix A is defined to be 

n 

det(A) = ^ sign(ai)sign(a2) 

fi,o-2eS'„ j=l 

We recall that biadjacency matrix A{x) = (aij) of a bipartite graph G = {V, W, E) is the \ V\ x \W\ matrix, 
defined as follows: We set 



aij :— 



if £ y- g ^ form an edge e with weight w(e), 
otherwise. 



1.2. Main results. 

Theorem 1. Let C he a linear binary code. Then there exists a tripartite triangular configuration A such 
that: // per(^A(a;)) — X^I^lo'^*^*' '^^^'"'^ A/\{x) is triadjacency matrix of A, then 

rn 
1=0 

where e is an integer linear in length of C . 

Proof. Follows from Theorems [51 [3] and 2] below. □ 

Theorem 2 (Rytif 9 ). LetC be a linear code over GF{p), where p is a prime. Then there exists a triangular 
configuration A such that: if X)i!lo '^^^^ weight enumerator of ker A then 

m 

Wc{x)=Y,a,x^^''''°^<''^'\ 

1=0 

where e is an integer linear in length of C. 

Theorem 3 (Ryti'f [5]). Let A 6e a triangular configuration. Then there exists a triangular configuration 
A' and weights w' : T(A') i— >■ M such that WkorA(2^) ~ PA'ix). 

Theorem 4. Let A be a triangular configuration with weights w : T(A) i— > M. Then there exists a tripartite 
triangular configuration A' and weights w' : r(A') K. such that Pa{x) — per(^A'(2:)) where A^i(x) is the 
triadjacency matrix of A' . 

Proof. Follows directly from Proposition 12.101 and Proposition 12.111 of Section [21 □ 

Definition 1.1. We say that an n x n x n 3-matrix A is Kasteleyn if there is 3- matrix A' obtained from A 
by changing signs of some entries so that per (A) = det(A'). 

Theorem 5. Let M be n x n matrix. Then one can construct m x m x m Kasteleyn 3-matrix A with 
m <n'^ + 2n and per(Af) — per(74). Moreover, Kasteleyn signing is trivial, i.e., per(j4) — det(A), and if M 
is non-negative then A is non-negative. 

Theorem [5l is proved in Section [3l 

2. Triangular configurations and permanents 

In this section we prove Theorem [H We use basic building blocks as in Rytff [B]. However, the use is 
novel and we need to stress the tripartitness of basic blocks. Hence we briefly describe them again. 

2.1. Triangular tunnel. Triangular tunnel is depicted in Figure [H An empty triangle is a set of three 
edges forming a boundary of a triangle. We call the empty triangles {a, b, c} and {a', b', c'} ending. 




Figure 1. Triangular tunnel 

Proposition 2.1. The triangular tunnel has exactly one matching with defect {a, 6, c} and exactly one 
matching with defect {a',6',c'}. □ 
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Figure 2. Tunnel tripartition 



Proposition 2.2. The triangular tunnel is tripartite. 
Proof. Follows from Figure [H 

□ 

2.2. Triangular configuration S^. Triangular configuration is depicted in Figured Letter "X" denotes 
empty triangles. We call these empty triangles ending. 




Figure 3. Triangular configuration 5*^ 

Proposition 2.3. Triangular configuration has one exactly perfect matching and exactly one matching 
with defect on edges of all empty triangles. 

Proof. The unique perfect matching is {ti, t2, ^4, is}- We denote it by M^{S^). The unique matching with 
defect on edges of all empty triangles is {is}. We denote it by M'^{S^). □ 

Proposition 2.4. Triangular configuration is tripartite. 
Proof. Follows from Figure HI 




Figure 4. Triangular configuration 5*^ with partitioning 

□ 
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Figure 5. Matching triangular triangle 



2.3. Matching triangular triangle. The matching triangular triangle is obtained from the triangular 
configuration and three triangular tunnels in the following way: Let Ti, T2 and T3 be triangular tunnels. 
Let , ; , q'^^ and t^^ ,r^^ be the ending empty triangles of Ti , T2 and T3 , respectively. Let if , i f ^ j ^ 
be ending empty triangles of S"^. We identify with if''; t'^^ with if and with tf''. The matching 
triangular triangle is defined to be riUS'^UT2UT3. The matching triangular triangle is depicted in Figure[5] 

Proposition 2.5. The matching triangular triangle has exactly one perfect matching and exactly 
one matching with defect {1, 2, 3, a, 6, c, a, /3, 7}. It has no matching with defect E, where 9 E C. 
{l,2,3,a,&, c,a,^,7}. 

Proof The perfect matching is M\S^)UM^{Ti)UM^{T2)UM^{T3). The matching is M°{S^)U 

M^(Ti) U M^{T2) U M^iTs). 

Any matching of the matching triangular triangle with defect E C {1, 2, 3, a, b, c, a, 13, -f} contains Mi(5^) 
or M'^{S'^). This determines remaining triangles in a matching with defect E C {1, 2, 3, a, &, c, a, /3, 7}. 
Hence, there are just two matchings and with defect E C {1, 2, 3, a, 6, c, a, /3, 7}. □ 

Proposition 2.6. Matching triangular triangle T is tripartite and there is a tripartition of T such that 
a,b,ce Ei; 1,2,3 € E2; a,i3,j e E3. 

Proof. Follows from Figure [4] and Figure [2] □ 

2.4. Linking three triangles by matching triangular triangle. Let A be a triangular configuration. 
Let ti, t2 and t^ be three edge disjoint triangles of A. 

The link by matching triangular triangle between ii,i2 and ^3 in A is the triangular configuration A' de- 
fined as follows. Let T be a matching triangular triangle defined in Section [^31 Let {a, b, c}, {1, 2, 3}, {a, /3, 7} 
be ending empty triangles of T. Let tl,tl, tf and tl,t2, i^ and ^3 , i§ , ^3 be edges of ti and t2 and ts , respectively. 
We relabel edges of T such that {a,b,c} = {t\,tl,tl} and {1,2,3} = {tl,tl,t^} and {a, ,3, 7} = 
We let A' := A U T. 

2.5. Construction. Let A be a triangular configuration and let w : T{A) 1— > R be weights of triangles. We 
construct a tripartite triangular configuration A' and weights w' : T(A') R in two steps. First step: We 
start with triangular configuration 

a; := Ai U A2 U A3 

where Ai, A2, A3 are disjoint copies of A. Let i be a triangle of A. We denote the corresponding copies of 
t in Ai, A2, A3 by ii,i2,i3, respectively. 

Second step: For every triangle t of A, we link ti,t2,i3 in by triangular matching triangle T. We 
denote this triangular matching triangle by Tt. Then we remove triangles ii,i2,i3 from A'l. We choose 
a triangle t' from M^{Tt) and set w'(t') := w{t). We set w'{t') := for t' G T{Tt) \ {t}. The resulting 
configuration is desired configuration A'. 

Proposition 2.7. Triangular configuration A' is tripartite. 



Proof. The triangular configuration A' is constructed from three disjoint triangular configurations Ai , A2 , A3 . 
From these configurations all triangles are removed. Hence, we can put edges E{Ai) to set Ei for i = 1,2, 3. 
The remainder of A' is formed by matching triangular triangles. Every matching triangular triangle connects 
edges of Ai, A2, A3. By Proposition 12 .61 the matching triangular triangle is tripartite and its ends belong to 
different partities. 

□ 

We recall 2^ denotes the set of all subsets X. We define a mapping / : 2-^^^^ H> 2-^''^ ^ as: Let 5 be a 
subset of T(A) then 

fiS) := {M\Tt)\t e 5} U {M\Tt)\t G T(A) \ S]. 

Proposition 2.8. The mapping f is a bijection between the set of perfect matchings of A and the set of 
perfect matchings of A' and w{M) — w'{f{M)) for every M C r(A). 

Proof. By definition, the mapping / is an injection. By Proposition 12. 5[ every inner edge of Tt, t G T{A), 
is covered by f{S) for any subset S of T{A). Let M be a perfect matching of A. We show that /(M) is 
perfect matching of A'. 

f{M) = {M\Tt)\t e M} U {M°{Tt)\t € T{A) \ M}. 

Let e be an edge of A and let ei, 62, 63 be corresponding copies in Ai, A2, A3. Let ti,t2, ■ . ■ ,ti be triangles 
incident with edge e in A. Let tk be the triangle from perfect matching M incident with e. By definition of 
A', the edges ei, 62, 63 are incident only with triangles of Tj. , i = 1, . . . ,1. The edges ei, 62, 63 are covered 
by M^{tk). The edges oi Tt-, i — 1, . . . ,1, i ^ k are covered by M^iTt^). Hence f{M) is a perfect matching 
of A'. 

Let M' be a perfect matching of A'. By Proposition M' = {M^{Tt)\t e 5} U {M°{Tt)\t £ T(A) \ S} 
for some set S. The set S* is a perfect matching of A. Thus, the mapping / is a bijection. 

□ 

Corollary 2.9. P/^{x) = Pa'{x). □ 

Proposition 2.10. Let A be a triangular configuration with weights w : T(A) n> R. Then there exist a 
tripartite triangular configuration A' and weights w' : T(A') H> R such that there is a bijection f between 
the set of perfect matchings 'P(A) and the set of perfect matchings ofV{A'). Moreover, w{M) — w'{f{M)) 
for every M G V(A), and Pa (a;) = P^'{x). 

Proof. Follows directly from Propositions 12.71 and 12.81 and Corollary 12.91 □ 

Proposition 2.11. Let A be a tripartite triangular configuration with tripartition Ei, E2, such that 
\Ei \ — |i?2| = I-E3I and let A/^{x) be its triadjacency matrix. Then Pa{x) — per(AA(a;)). 

Proof. We have 

n 

fl ,'^2 GS„ i=l 

Every perfect triangular matching between partities Ei, E2, E3 can be encoded by two permutations cti,(72 
and vice versa. If matching M is a subset of T(A), then 

n n 
i=l i=l 

where [ijk] denotes a triangle of A with edges i,j, k. If M is not a subset of T(A), then there is i such that 
ata,{i)a2{r) = 0. Hcncc J]^"^^a;"'([^'^l(0'^2(^)]) = 0. Therefore 

n 

□ 



6 



3. Kasteleyn 3-matrices 



We first introduce a necessary condition for a 3-matrix to be Kasteleyn. Let A be a |Vo| x \Vi\ x IV2I 
non-negative 3-matrix, where \Vi\ = m,i = 1,2,3. We first define two bipartite graphs Gi,G2 as follows. 
We let, for i = 1,2, Gf = Gi = (Vb, Vi, Ei) where 

El = {{a, b}\a G Vo,b £ Vi and Aabc for some c} , 

and 

E2 = {{a, c}\a e Vo, c S V2 and Aabc for some b} . 

Theorem 6. If A is such that both are Pfaffian bipartite graphs then A is Kasteleyn. 

Proof. Let Mi be the biadjacency matrix of Gi and let sign^ : E{Gf) ^ {—1, 1} be the signing of the entries 
of Mi which defines matrix such that per(Mj) = det(Mj'). We define 3-matrix A' by 

Kbc = signi({a,6})sign2({a,c})^a6c- 

We have 

det(^') =X!^^Sn(CTi) X ^sign((72) ]Jsign2({i,f72(i)})[signi({j,ai(j)})AjVi(j)<r20)]- 

By the construction of siguj we have that for each (T2 and each cti, if Y[j ^jcri(j)cr2(j) 7^ then 

sign((72)JJ signs ({j, (72 (j)}) = 1. 
j 

Hence 

det(A') = ^sign(CTi) X ^J|[signi({j,cri(j)})Aj>i(j)^,(j)] 

o"l 0-2 j 

= XI ^ 51 sign(o-i) n^ig^i 0')})^J>iO)<^20)- 

(T2 <Tl j 

Analogously by the construction of siguj^ we have that for each ai and each a2, if Yij ^jcriU)c2U) then 

sign((7i) J|signi({j,(7i(j)}) = 1. 
j 

Hence 

det(A') = X =per(A). 

□ 

In the introduction we defined the triadjacency 3-matrix of a triangular configuration as the adjacency 
matrix of the edges of the triangles. We also defined a matching of a triangular configuration as a set of 
edge-disjoint triangles. In this section it is advantageous to consider any 3-matrix with 0, 1 entries as the 
adjacency matrix of vertices of a triangular configuration. Hence we need the following notions. 

A triangular configuration A is vertex-tripartite if vertices of A can be divided into three disjoint sets 
yi,V2, V3 such that every triangle of A contains one vertex from each set Vi, V2, V3. We call the sets Vi, V2, V3 
vertex-tripartition of A. 

The vertex-adjacency 3-matrix A{x) = (aijk) of a tripartite triangular configuration A with vertex- 
tripartition Vi,V2,l^ is the 1 X |\^| x IV3I three dimensional array of numbers, defined as follows: We 

set 

{j.w(t) jf j g j g /j g ^ forms a triangle t with weight w(t), 
otherwise. 

We will need the following modification of the notion of a matching. A set of triangles of a triangular 
configuration is called strong matching if its triangles are mutually vertex-disjoint. 
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Proof of Theorem\^ Let M he & n x n matrix and let G = {Vi, V2,E) be the adjacency bipartite graph of 
its non-zero entries. We have \Vi\ = IV2I = n. We order vertices of each Vi, i — 1,2 arbitrarily and let 
Vi — {v{i, 1), . . . , v{i, n)}. Let V- = {v'{i, 1), . . . , v'{i, n)} be disjoint copy of Vi, i = 1,2. 

We next define three sets of vertices Wi, W2, Wq and system of triangles A(G) = A so that each triangle 
intersects each Wi in exactly one vertex. 

Wi =ViUVlu{w{l,e)\ee E}, 

W2^V2LI K,' U {w{2, e)\e G E}, 

Wo = {w{0,e)\e e E} U {wiO,i,j)\i = 1,2; j ^l,...,n}, 

A = Ue=abeE {{a, b, w(0, e)), (w(0, e), e), w(2, e))}U 

Uj^i {{wi0,l,j),v'i2,j),w{l,e))\v{l,j) e e}U 

Uj^i {{w{0,2,j),v'{l,j),w{2,e))\v{2,j) G e}. 

We let A be the vertex-adjacency 3-matrix of the triangular configuration T(G) = T ~ {Wq, Wi, W2, A). 
We first observe that both bipartite graphs Gi, G2 of this triangular configuration are planar; let us consider 
only Gi, the reasoning for G2 is the same. First, vertices v'{l,j) and w{0,2,j) are connected only among 
themselves in Gi. Further, the component of Gi containing vertex v{l,j) contains also vertex w(0, l,j) and 
consists of degG{v{l,j)) disjoint paths of length 3 between these two vertices. Here degG{v{l,j)) denotes 
the degree of v{l,j) in graph G, i.e., the number of edges of G incident with v{l,j). Thus, by Theorem[Bl 
A is Kasteleyn. 

We next observe that Kasteleyn signing is trivial. Let Di be the orientation of Gi in which each edge is 
directed from Wq to Wi . In each planar drawing of Gi , each inner face has an odd number of edges directed 
in Di clockwise. This means that Di is a Pfaffian orientation of Gi, and per(A) — det{A) (see e.g. Loebl [5] 
for basic facts on Pfaffian orientations and Pfaffian signings). 

Finally there is a bijection between the perfect matchings of G and the perfect strong matchings of 7": if 
P C i? is a perfect matching of G then let 

P(r) = {(a,6,w(0,e))|e = a6eP}. 

We observe that P{T) can be uniquely extended to a perfect strong matching of T, namely by the set of 
triples S*! U 5*2 U 6*3 where 

Si = UegB\p{(u;(0, e), w{l, e), w{2, e))}, 

^2 - \J';^A{w{0,l,j),v'{2,j),w{l,e))\v{l,j) e e e P}, 

^3 = U;U{(i«(0,2,i),i;'(l,j),u;(2,e))|i;(2,j) e e e P}. 

Set S'l is inevitable in any perfect strong matching containing P{T) since the vertices w{Q, e); e ^ P must 
be covered. This immediately implies that sets 5*2 , 6*3 are inevitable as well. 

On the other hand, if Q is a perfect strong matching of T then Q contains P(7~) for some perfect matching 
P of G. 

□ 

4. Application to 3D dimer problem 

Let Q be cubic n x n x n lattice. The dimer partition function of Q, which is equal to the generating 
function of the perfect matchings of Q, can be by Theorem [5] identified with the permanent of the Kasteleyn 
vertex- adjacency matrix of triangular configuration T{Q). Natural question arises whether this observation 
can be used to study the 3D dimer problem. 

We first observe that the natural embedding of Q in 3-space can be simply modified to yield an embedding 
of T{Q) in 3-space. This can perhaps best be understood by figures, see Figure El this figure depicts 
configuration T{Q) around vertex v oi Q with neighbors ui, . . . ,uq. 

Triangular configuration T{Q) is obtained by identification of vertices t;^, i = 1, . . . , 6 in the left and right 
parts of Figure [H Now assume that the embedding of left part of Figure [6] is such that for each vertex v of 
Q, the vertices vi, . . . ,vq belong to the same plane and the convex closure vi, . . . ,vq intersects the rest of 
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Figure 6. Configuration T{Q) around vertex u of Q with neighbors mi, . . . , ue so that «i, U3, ^4, ue belong 
to the same plane in the 3-space, U2 is 'behind' this plane and ur, is 'in front of this plane. Empty vertices 
belong to Wo, square vertices belong to Wi and full vertices belong to W2- 

the configuration only in wi, . . . , wg. Then we add the embedding of the right part, for each vertex v of Q, 
so that xi belongs to the plane of the ViS and X2 is very near to xi but outside of this plane. 

Summarizing, the dimer partition function of a finite 3-dimensional cubic lattice Q may be written as the 
determinant of the vertex-adjacency 3-matrix of triangular configuration T{Q) which preserves the natural 
embedding of the cubic lattice. Calculating the determinant of a 3-matrix is hard, but perhaps formulas for 
the determinant of the particular vertex-adjacency 3-matrix of T{Q), illuminating the 3-dimensional dimer 
problem, may be found. An example of a formula valid for the determinant of a 3-matrix is shown in the 
next subsection. It is new as far as we know but its proof is basically identical to the proof of Lemma 3.3 of 
Barvinok [T]. 

4.1. Binet-Cauchy formula for the determinant of 3- matrices. We recall from the introduction that 
the permanent of a n x n x n 3-matrix A is defined to be 

n 

per(A) = ^ n^^'^iW'^^W- 

0'l,0-2eS„ 4=1 

The determinant of a n x n x n 3-matrix A is defined to be 

n 

det(A) = ^ sign(cri)sign(cr2)n 

0'1,0'2GS„ 4=1 

The next formula is a generalization of Binet-Cauchy formula (see the proof of Lemma 3.3 in Barvinok |T|). 

Lemma 4.1. Let A^, A^, A'^ be real r x n matrices, r < n. For a subset I C {1, . . . ,n} of cardinality r we 
denote by Aj the r x r submatrix 0/ the matrix A'^ consisting of the columns of indexed by the elements 
of the set I. Let C be the 3-matrix defined, for all ii, 12, is by 

n 

J=l 

9 



Then 

det(C) = ^per(A})det(A?)dct(yl?), 
/ 

where the sum is over all subsets / C {1, . . . , n} of cardinality r 
Proof. 

r n 

det(C) = J2 sign(ai)sign(a2) J] E ^hK(^),J^U^hJ 

r 

= J2 sign(cri)sign((72) X ^ H ^L.^^iW.j.^^sW 

r 

Now, for all J — (ji, . . . , j^) we have 

r 

J2 sign(cri)sign(cr2) x H ^L-.^^i(«),i.^^2(«) j. 

=(n 4.j(E «ig^(^i) n <w,.j(E «ig^(^2) n 

i — 1 (Ti (72 2—1 

r 

=(n4.Jdet(i2)det(i3), 

where Aj denotes the r x r matrix whose ith column is the j'^th cohimn of matrix . 

If sequence J contains a pair of equal numbers then the corresponding summand is zero, since det(Aj) is 
zero. Moreover, if J is a permutation, and J' is obtained from J by a transposition, then 

det(i2 ) det(ii}) = det(i^,)det(ii},)- 

Therefore Lemma 14^1 follows. 

□ 

References 

[1] A. I. Barvinok. Now algorithms for linear k-matroid intersection and matroid k-parity problems. Mathematical Program- 
ming, 69:449-470, 1995. 

[2] P. W. Kasteleyn. The statistics of dimers on a lattice. I. The number of dimer arrangements on a quadratic lattice. Physica, 
27(12):1209-1225, 1961. 

[3] P. W. Kasteleyn. Dimer Statistics and Phase Transitions. Journal of Mathematical Physics, 4(2):287-293, 1963. 

[4] P. W. Kasteleyn. Graph theory and crystal physics. In Graph Theory and Theoretical Physics (ed. by F. Harary), Academic 

Press, New York, pages 43—110, 1967. 
[5] M. Loebl. Discrete Mathematics in Statistical Physics. Vieweg+Teubner, 2010. 
[6] G. Polya. Aufgabe 424. Arch. Math. Phys. Ser., 20:271, 1913. 

[7] N. Robertson, P. D. Seymour, and R. Thomas. Permanents, Pfafhan orientations, and even directed circuits. Math. Ann., 
150:929-975, 1999. 

[8] P. Ryti'i'. Geometric representations of binary codes and computation of weight enumerators. Adv. in Appl. Math., 45:290- 
301, 2010. 

[9] P. Rytff. Geometric representations of linear codes, submitted, 2011. 
[10] H. N. V. Temperley and M. E. Fisher. Dimer problem in statistical mechanics-an exact result. Phil. Mag., 6:1061—1063, 
1961. 



10 



